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Recent theoretical work has shown that so-called pulled fronts propagating into an unstable 
state always converge very slowly to their asymptotic speed and shape. In the the light of these 
predictions, we reanalyze earlier experiments by Fineberg and Steinberg on front propagation in a 
Rayleigh-Benard cell. In contrast to the original interpretation, we argue that in the experiments the 
observed front velocities were some 15% below the asymptotic front speed and that this is in rough 
agreement with the predicted slow relaxation of the front speed for the time scales probed in the 
experiments. We also discuss the possible origin of the unusually large variation of the wavelength 
of the pattern generated by the front as a function of the dimensionless control parameter. 



I. INTRODUCTION 

Although the propagation of a front into an unstable 
state plays an important role in various physical situ- 
ations ranging from the pearling instability (l],|^] to di- 
electric breakdown j| , detailed experimental tests of the 
explicit theoretical predications, especially those for the 
velocity of so-called "pulled" fronts are scarce. One of 
the reasons lies in the difficulty in preparing the system 
in the unstable state. 

If the initial front profile is steep enough the prop- 
agating front converges to a unique shape and veloc- 
ity. Theoretically, one distinguishes two regimes for front 
propagation into unstable states: the so-called "pushed" 
regime, where the front is driven by the nonlinearities 
and the so-called "pulled" regime where the asymptotic 
velocity of the propagating front, v as , equals the spread- 
ing speed v* of linear perturbations around the unstable 
state: v as — v*. Pushed fronts are by definition those 
for which the asymptotic speed w as is larger than v*: 
v as > v* . It is thus as if a "pulled" front is literally 
"pulled" by the leading edge whose dynamics is driven 
by linear instability of the unstable state [^j-Q ; the non- 
linearities merely cause saturation behind the front. We 
focus here on the experimental tests of the dynamics of 
such pulled fronts; since v* is determined by the equa- 
tions linearized about the unstable state, the front ve- 
locity of pulled fronts can often be calculated explicitly, 
even for relatively complicated situations. 

There have been two experiments aimed at testing the 
predictions for the speed of pulled fronts. Almost 20 
years ago Alders and Cannell || studied the propagation 
of a vortex- front into the laminar state in rotating Taylor- 
Couette flow. The measured velocities were about 40% 
smaller than expected from the theoretical predictions. A 



few years later, however, Fineberg and Steinberg (FS) jjj] 
published data which appeared to confirm the expected 
velocity in a Rayleigh-Benard convection experiment to 
within about 1%. The issue then seemed to be settled 
when it was also shown that the discrepancy observed by 
Ahlers could be traced back to slow transients 0] . 

The theoretical developments of the last few years give 
every reason to reconsider the old experiments by FS: It 
has been shown that the convergence of the velocity 
of pulled fronts is always very slow, in fact with leading 
and subleading universal terms of 0(l/t) and 0(l/t 3 / 2 ) 
with prefactors which follow from the linearized equa- 
tion. This slow relaxation implies that it will in general 
be very difficult to measure the asymptotic front speed to 
within a percent or so in any realistic experiment. Hence, 
from this new perspective, the proper question is not why 
in the Taylor-Couette experiment the measured velocity 
was too low, but why apparently in the Rayleigh-Benard 
experiment of FS the asymptotic front speed was mea- 
sured. 

The main purpose of this paper is to address this is- 
sue, and to reanalyze the experiments in the light of the 
present theory. We will conclude that the data of FS 
actually do show signs of the predicted power law con- 
vergence of the front velocity to an extrapolated asymp- 
totic value which is about 15% larger than their tran- 
sient value. This of course implies that there is then a 
discrepancy of order 15% between the value of v as = 
as extrapolated from their data, and the one claimed in 
the original experiments. We will argue that the most 
likely reconciliation of the two results is that the value of 
the correlation length £o in the experimental cell of FS 
is somewhat larger than the theoretical value used by FS 
to interpret their data. 

Of course, only new experiments can settle whether the 
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interpretation we propose is the correct one. We do con- 
sider new experiments along the lines of FS in fact very 
desirable, not so much as they might settle the numeri- 
cal value of the velocity, but more because they hold the 
promise of being the first accurate experimental test of 
the universal power law relaxation of pulled fronts. 

In section II we will first summarize the relevant theo- 
retical predictions for the velocity of pulled fronts. Then 
we will discuss the experiments of FS in the light of these 
results in section III, where we will also reanalyze their 
data. Finally, in section IV we turn to a brief discussion 
of the wavenumber of the pattern selected by the front. 
Here, the results of FS were not quite consistent with 
the predictions for the asymptotic wavenumber from the 
Swift- Hohcnberg equation. As we shall discuss, the wave- 
length of the pattern is affected by various effects which 
are not easily controlled, but the most likely interpreta- 
tion of the data of FS is that they did not observe the 
asymptotic wavelength behind the front, but the local 
wavenumber in the leading edge of the front. Indeed it 
is in general difficult to test the theory by studying the 
asymptotic pattern wavelength and the convergence to 
the asymptotic value experimentally. 



II. SUMMARY OF THEORETICAL 
PREDICTIONS 

A. Asymptotic speed and power law convergence 

Just above the onset of a transition to stationary fi- 
nite wavelength patterns, for small dimensionless control 
parameters e the slow dynamics on length scales larger 
than the wavelength of the pattern can be described by 
the Ginzburg-Landau amplitude equation jll]|l^] 



T d t A = eA + $d 2 x A-g\A\ 2 A 



(1) 



The time scale To and length scale £o a s well as the non- 
linear saturation parameter g depend on the particular 
system under study. 

The asymptotic spreading speed v* of linear pertur- 
bations around the unstable state is in general obtained 
from the linear dispersion relation ui(k) of a Fourier mode 
e -iut+ikx through 
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This yields for the Ginzburg-Landau equation [[L3| 
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In the experiments with which we will compare, the 
scaled velocity 
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is often used. According to (|3|), for a pulled front the 
asymptotic value v as — 2. 

The above results were known in the eighties, at the 
time when the experiments were done. The crucial in- 
sight of the last few years is the finding that the con- 
vergence or relaxation towards the asymptotic velocity 
v* of pulled fronts is always extremely slow: the general 
expression for the time dependent velocity v(t) emerging 
from steep initial conditions (i.e., decaying faster than 
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For the case of the Ginzburg-Landau amplitude equation, 
we then have 
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It is important to stress that the above expression for 
v(t) is exact but asymptotic — this is illustrated by the 
fact that at time et/ro = n the subdominant i~ 3 / 2 term 
is equal to the first correction term of order i -1 in abso- 
lute value, but of opposite sign. Thus, although for any 
sufficiently long time i, the above expression will always 
become accurate, at any finite time, however, the expres- 
sion might only yield a good estimate. In fact, in practice 
one usually has to go to dimensionless times et/ro of or- 
der 10 or larger for the asymptotic expression to become 
accurate, while for dimensionless times et/ro m the range 
3 to 10 the first correction term yields a reasonable order- 
of- magnitude estimate [Q. As we shall see below, in the 
Rayleigh-Benard experiments || the maximum dimen- 
sionless time et/ro that can be probed is about 3 to 4. 
In comparing with experiments and in making order of 
magnitude estimates, we will therefore only use the first 
correction term. Finally it is important to realize that 
after how long a time these expressions become really 
accurate, depends also on the initial conditions. 



B. Dependence on initial conditions 

In order to illustrate how accurate these expres- 
sions are in practice, we numerically integrate the real 
Ginzburg-Landau equation starting with an exponen- 
tially decaying initial condition with steepness /i [13 
A(x,t = 0) = ae~^ x . The result is shown in Fig 
for various values of \i. We see that initially the veloc- 
ity falls off quickly and then approaches the asymptotic 
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velocity from below, and that the asymptotic expression 
(||) in practice does yield a reasonable estimate of the 
time dependent velocity for values of the scaled time of 
order 3 and larger. 

We also note that the theoretical analysis shows that 
for initial profiles falling off exponentially with steepness 
fi < ji* , the asymptotic velocity lies above v* and is given 
by v as (fi) = fi + — . The dotted line in Fig. |l| shows an 
example of a case with fi slightly less than ji* , for which 
v as = 2.05. In this case, the time-dependent scaled ve- 
locity is approximately equal to 2 at times of order 3 to 
4. 

In the experiments, the typical protocol was to simul- 
taneously increase the heat flux from an initial state at 
£i < (with a typical value of e$ = —0.015) to a su- 
percritical value e/, and switching on the end heater. A 
rapid convergence to the asymptotic velocity due to spe- 
cial initial conditions with A ~ A* would have required 
setting 6i ~ —ef for all e/ and would probably also have 
required switching on the end heater before the value of 
e was changed, so that effectively at the end a convection 
pattern was prepared whose amplitude decayed exponen- 
tially into the bulk. Thus, it is theoretically possible 
to select special initial conditions so as to get a scaled 
velocity around 2 at some finite time, but the finetun- 
ing necessary to do so is so sensitive that we consider it 
very unlikely that experimental observations done over a 
range of values of e are due to initial condition effects. 
As stressed before, only new experiments can completely 
rule out this possibility, however. 
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FIG. 1. Velocity of fronts propgating propgating into 
the unstable state A = of the real Ginzburg-Landau 
equation with exponentially decaying initial conditions 
A = \J e /9 e ~^ x ■ The velocity is obtained by tracing the 
position of the point where A reaches half its asymptotic 
value. Since e and £o just set the length and time scale, 
v(t) plotted as a function of the scaled time et/ro is pa- 
rameter-independent. Different initial conditions or tracing 
a different value of A to determine the velocity yield different 
transient behavior, but for sufficiently steep initial conditions, 
all curves converge to the analytic formula for late times. 



III. REEXAMINATION OF THE 
RAYLEIGH-BENARD FRONTS 



In the long quasi-one-dimensional Rayleigh-Benard cell 
of FS H the front propagation was initiated by simulta- 
neously increasing the heat flux to a supercritical value 
e > (e = -^f- — 1) and switching on a heater at the 
ends of the long cell. A vortex front is induced near the 
this heated end-wall and the propagation of this front 
into the unstable conductive state is then studied. Thin 
fins were attached to the long sides of the cell to avoid 
both induction of long rolls and pinning. Both because 
of the fact that the initial perturbation was caused by 
heating at the end and because the state before bringing 
the temperature difference beyond its critical value was 
unrelated to the final value of e, there is every reason 
to believe that the experimental protocol did not create 
any special initial conditions that can not be considered 
sufficiently steep or localized. 

As FS point out, since the front velocity v* grows as 
y/e while the growth rate in the bulk grows as e, fronts 
can only be observed and in fact dominate the dynamics 
for small enough e. In practice the pattern could be dis- 
tinguished from the bulk noise up to a time t\, g — riToe 
where the numerical factor n is of order 3-4. This deter- 
mines some upper limit e on e f° r which the front can 
advance of the order of (l/m)th of the cell length I before 
bulk growth takes over: 

Let us now estimate, using our analytical estimate (^), 
the relative importance of the correction term at the lat- 
est times of order t^ g — riT^e -1 at which measurements 
can be taken. Substitution gives 

«(^)« 2 (!-^)' ( 10 ) 

Thus, for the latest time accessible in the experiment, we 
obtain with the empirical experimental value n =3-4 a 
velocity which is of order 20 ± 5% below the asymptotic 
value. Although the asymptotic formula may not be ac- 
curate yet at such early times, the numerical results of 
Fig. [I] lead one to expect corrections of the same order of 
magnitude: these times also correspond to scaled times 
of order n in the numerical simulation plots of Fig. [l], 
and as we have seen, over this time range, the velocity 
is also suppressed by about 15% relative to the asymp- 
totic value. Hence, by any reasonable estimate, the slow 
relaxation can not be negligible in the experiments! 

As stated in the introduction, FS quoted a measured 
velocity v = 2.01 ± 0.02 for their experiments, contrary 
to our expectation that their scaled velocity should have 
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been in the neighborhood of 1.7. According to the the- 
oretical results, one should expect to get to within 1% 
of the asymptotic velocity only around a dimensionless 
time of order 100! 

In our view, the most plausible explanation for the ori- 
gin of this discrepancy is that the value £o was actually 
larger in the experiments than the value used in the anal- 
ysis. The value of tq was experimentally confirmed to be 
very close to the theoretical value |l4|,[l5|]. For £q, how- 
ever, the theoretical value was used without independent 
experimental check [Q . Because of the special design of 
the cell with side-fins to create one-dimensional patterns, 
a different value of £o might not be unexpected. In fact, 
an indication that £o m the experiments was larger than 
the theoretical value used in the analysis, comes from the 
observation of FS that the value of the wavelength at on- 
set was 13% larger than the theoretical value. This might 
indicate that all lengths in the experiments are a factor 
1.13 larger than the theoretical values, and this is pre- 
cisely the factor needed to reconcile the front data with 
the theoretical expectations! One should keep in mind, 
though, that £ is determined by the the curvature of 
e versus k tongue around the critical wavenumber, and 
that it is not guaranteed that both are changed by the 
same factor; only independent measurements can fully 
settle this issue. 

We now show that a reanalysis of the data of FS actu- 
ally gives quite convincing evidence for slow convergence 
effects in the experimental fronts. 

FS measured the velocity by comparing the front with 
itself at various time intervals, by appropriately shifting 
the traces back. This yielded a set of points in the Aa:, At 
plane which appeared to lie on a straight line. However 
a possible relaxation of the velocity was masked since 
points from early and later times will approximately fall 
in the same place in the plane, and because the front 
shape also has an asymptotic 1/t relaxation. We there- 
fore have tried to reanalyze the raw data of figure 2 from 
the FS paper; this space-time plot of a propagating from 
is reproduced in the top panel of our Fig. ^. We define 
the position of the front as the point where the interpo- 
lation of the maxima of the profile equals some fraction 
of its maximum in the bulk (we chose 0.4). Our data for 
v(t) obtained this way are shown in the lower panel of 
Fig. ^. Whereas the local velocity initially slightly de- 
creases, an increase for dimensionless times larger than 
0.5 is evident. We stress that this qualitative behavior 
is independent of the choice of the parameters. In order 
to compare quantitatively to the predictions for the re- 
laxation, we have used the value for tq given by FS but 
increased the value of £o by 13% on the basis of the ar- 
gument given above. Clearly, with this choice, the data 
are certainly consistent with the analytical as well as nu- 
merical estimates of the velocity relaxation — in fact, in 
a way the data are the first experimental indications for 
the universal power law relaxation of pulled fronts U7j . 
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FIG. 2. Top panel: shadowgraph trace of a propagating 
front in the experiments of FS for e = 0.012 H. The time 
difference between successive traces is 0.42 t v , where t v is the 
vertical diffusion time in the experiments, and the distances 
are measured in units d (the cell height), (from Middle 
panel: similar data obtained from numerical integration of 
the Swift-Hohenberg equation also at e = 0.012 starting with 
a localized initial condition. The time difference between suc- 
cessive traces corresonds to 0.42 t v . Bottom panel: Velocity 
versus time in the experiment, as obtained by interpolating 
the maxima of the traces in the top panel, as explained in the 
text. The dashed line shows the analytical result (Q) and the 
dotted curve the result of the amplitude equation simulation 
of Fig. [l] with fi/fj.* — 1.2. Note that the curves are not fitted, 
only the absolute scale is affected by adjusting £o- 
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IV. RELAXATION OF WAVELENGTH 



FS also studied the problem of the selection of the 
wavelength A. As we discussed before, the actual value 
of the wavelength A c of the patterns is at criticality about 
13% off from the theoretical value; however, we are not 
interested here in the absolute value, but in the relative 
variation of A c /A. 

The difficulty of comparing theory and experiment on 
the variation of the wavelength is that the only theor- 
ically sharply defined quantity is the wavelength suffi- 
ciently far behind the front, A as , and that one has to 
go beyond the lowest order Ginzburg-Landau treatment 
to be able to study the pattern wavelength left behind. 
E.g., if we use a Swift- Hohcnbcrg equation for a system 
with critical wavenumber k c and bare correlation length 

d t u= —{1 + ^) u + eu-u, (11) 




e 

FIG. 3. The selected wave number A c /A as function of e. 
The data points of FS are denoted by circles, our numeri- 
cal results on the Swift-Hohenberg by triangles. The dashed 
lines shows the prediction for the asymptotic wavelength, the 
dashed-dotted line shows the analytic result with relaxation 
terms at t = 6ro/^/e and the grey line shows the result 
for the local wavelength in the leading edge. 



then a node counting argument [QJ^] yields for the asymp- 
totic wavelength A as far behind the front || : 
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In the Rayleigh-Benard experiments, k c « 2.75/ d, where 
d is the cell height; the theoretical value is £o = 0.385d, 
so if our conjecture that the value is some 15% larger is 
correct, we get £o ~ 0.44d. This then gives 

w 1 + 0.34e. (13) 

Aas 

As we stressed already above A os is the wavelength far 
behind the front; for a propagating pulled front, there is 
another important quantity which one can calculate ana- 
lytically, the local wavelength A* measured in the leading 
edge of the front. For the Swift-Hohenberg equation, one 
gets for this quantity |^,|6| 

~ 1 + ^e (e<<1) 

wl + e, (14) 

where in the last line we have used the experimental val- 
ues. 



Let us now discuss the experimental findings in the 
light of these results. In Fig. || we show both the exper- 
imental values of X c /X eX p as determined experimentally 
by FS, and results for the Swift-Hohenberg equation with 
the value (|l3|) relevant for the experiments. 

(i) For small e, the experimental values are roughly 
linear in e, but when a fit is made over the whole range 
of e values studied, a square root behavior, as proposed 
by FS, would probably be better. 

(ii) The experimental values for the wavelength ratio 
deviate about a factor of 3 from the theoretically ex- 
pected value for the ratio far behind the front, A c /A as . 

(Hi) Just like the front speed converges very slowly 
to the asymptotic value, so does the local value of the 
wavelength behind the front The relaxation of the 
velocity is plotted with a dashed-dotted line in Fig. || for 
times about 6to / y/e, which is the time it takes for a front 
to propagate to close to the center of a system about as 
large as the experimental cell. Clearly, the wavelength ra- 
tio due to slow relaxation lies below the asymptotic value, 
and hence further away from the experimental data for 
small e. Also, numerical results for the wavelength of the 
fourth "roll" measured at the same time (indicated by 
the dotted line, lie below the asymptotic curve. 

(iv) FS have measured the wavelength for very small 
values of e, down to 4 1CP 4 . However, already for values 
as small as 0.01, the coherence length £ = £o/i/e is about 
2A C in the experiments. The total front width is several 
times this number, and for even smaller values of e the 
front width is even smaller. Since the total experimental 
cell was about 12A C long in the experiments, it is clear 
that over much of the small e range, one would not ex- 
pect to see well-developed fronts in the experiments. In 
other words, it is very unlikely that over the experimen- 
tal range of e-values, one has a chance to measure A as of 
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a well-developed pattern behind a front. 

(v) The grey line in the plot shows our analytical re- 
sult for the wavelength ratio in the leading edge of the 
front. Clearly, this line follows the data for small e quite 
well. In view also of point (iv) above that it will be hard 
to obtain well-developed fronts for small e, we propose 
as a tentative explanation of the data that in the small 
e range, one actually measures the emergent roll pattern 
associated with the leading edge of a front. Of course, 
only new experiments can decide on the validity of this 
suggestion. 

(vi) We mention that the variation of the wavelength 
ratio with e depends also on the third order derivative 
term in the expansion of the dispersion relation around 
the critical wavenumber. This term is not modeled cor- 
rectly in the Swift-Hohenberg equation, but may have to 
be taken into account in a full comparison of theory and 
experiment. 

(vii) We finally mention that in the experiments there 
was an up-down asymmetry in the rolls. We have inves- 
tigated whether this could be a source of the discrepancy 
between the asymptotic wavelength ratio and the ob- 
served one, by studying a Swift-Hohenberg equation with 
a symmetry-breaking quadratic term. However, with this 
term, the wavelength ratio apears to decrease away from 
the experimental values. 



V. CONCLUSION 

It was recently discovered that quite generally pulled 
fronts relax very slowly to their asymptotic velocity. 
Comparison of the experimental data for the velocity 
with numerical simulations and analytical estimates give, 
in our view, evidence that these experiments provide 
clear signs of the presence of such slow relaxation effects, 
although the time scales that can be probed experimen- 
tally are too short to test the general power law relax- 
ation. Theoretically the only other viable option to rec- 
oncile theory with the interpretation originally proposed 
H is that somehow special initial conditions created an 
initial convection profile with precisely the right spatial 
decay into the bulk. As we discussed in section III, in 
our view this is an unlikely interpretation, but only new 
experiments can settle this issue completely. 

While measurements of the wavelength of the pattern 
generated by a front are even more difficult to interpret 
than those of the velocity, our analysis indicates that 
in the small-e regimes a well-developed front does not 
fit into the experimental cell, and that as a result one 
probes the local wavelength in the leading edge of the 
front rather than the well-developed asymptotic wave- 
length behind it. The analytical estimates are consistent 
with this suggestion. 

We hope that this work will trigger new experimental 
activity to investigate these issues — experiments along 
these lines hold the promise of being the first ones to see 



the universal power law relaxation of pulled fronts. 
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